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OBJECTIVE 

Gather  data  on  magnetospheric  substonns,  and  the  ensuing  ionospheric  storms  with 
» a view  to  the  development  of  models  useful  in  minimizing  the  detrimental  impact  of  these 

processes  on  Naval  communications. 

« 

RESULTS 

1 . The  general  morphology  of  both  magnetospheric  substorms  and  ionospheric 
storms  has  been  documented,  although  much  more  study  is  necessary  before  any  accurate 
modeling  can  be  done. 

2.  A model  of  the  “shock  front"  preceding  some  ionospheric  storms  has  been 
developed. 


RECOMMENDATIONS 

Additional  data,  both  magnetospheric  and  ionospheric,  must  be  acquired  before 
accurate  and  detailed  models  of  these  processes  can  be  developed.  Presently,  a large-scale 
effort  is  being  pointed  this  direction,  in  the  International  Magnetospheric  Study  (IMS). 
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PREFACE 


Geomagnetic  and  ionospheric  variations  are  extremely  important  to  the  Navy  because 
many  and  surveillance  systems  depend  upon  the  stability  of  the  geomagnetic  field  and 
ionosphere  for  reliable  operation.  For  example,  rapid  geomagnetic  pulsations  may  obscure 
submarines  from  our  ASW  detectors,  and  ionospheric  storm  depletions  or  enhancements  of 
the  F-region  may  compromise  covert  Naval  communications.  To  alleviate  the  damaging 
effects  of  these  variations,  a real-time  environmental  prediction/assessment  system  must  be 
developed  with  inputs  into  a model  describing  effects  of  the  variations  on  the  system. 

This  report  documents  the  present  understanding  of  both  magnetospheric  substorms 
and  ionospheric  storms,  and  suggests  a program  to  expedite  the  development  of  an  opera- 
tional model  for  Naval  use. 

Magnetospheric  substorms  appear  to  be  the  causative  mechanism  of  the  global  iono- 
spheric storm.  A digest  of  the  major  models  of  magnetospheric  substorms  is  presented  in 
chapter  1.  Chapter  2 discusses  the  ionospheric  storms  that  are  identified  with  magneto- 
spheric substorms. 
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CHAPTER  1 : MAGNETOSPHERIC  SUBSTORMS 
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INTRODUCTION 

The  solar  wind  (a  neutral  plasma)  flowing  by  and  interacting  with  the  earth’s  mag- 
netic field  undergoes  charge  separation.  Thus,  the  mechanical  energy  of  the  solar  wind  is 
converted  to  an  electric  potential  by  the  geomagnetic  field  in  a dynamo-like  manner. 
Akasolu  (1975)*  estimates  that  the  solar  wind  is  putting  energy  into  the  magnetosphere  at 
the  rate  of  ~10'tJ  erg/s.**  Assuming  that  the  magnetosphere  cannot  hold  an  infinite 
amount  of  energy,  we  might  expect  some  mechanism  for  energy  release.  A magnetospheric 
substorm  is  an  implosive  release  of  this  magnetospherically  stored  energy  into  the  polar 
ionosphere,  and  can  be  viewed  as  a form  of  capacitive  discharge  through  a current  system 
to  be  described  later. 

The  energy  release  (or  conversion  of  potential  energy  into  mechanical  energy)  is 
manifested  through  many  forms  in  the  polar  ionosphere,  as  documented  by  Akasolu 
(1968):*** 


Auroral  substorm.  Greatly  enhanced  auroral  activity  - specific  patterns  of 

variation 

Polar  magnetic  substorm.  Enhanced  auroral  electrojet,  creating  magnetic 

disturbances 

Ionospheric  substorm.  Anomalous  ionization  and  heating  in  lower  ionosphere 
by  auroral  particles,  generating  traveling  wave  disturbances  and  infrasonic  waves  — resulting 
in  major  changes  in  circulation 

X-ray  substorm.  Bremsstrahlung  from  energetic  precipitated  electrons 

Proton  aurora  substorm.  Precipitating  low-energy  protons  enhance  auroral 
luminosity  in  the  Balmer  lines 

Vlf  emission  substorm.  Hiss  and  chorus  coming  from  precipitated  particles 

Micropulsation  substorm.  Ulf  electromagnetic  waves,  thought  to  be  of  mag- 
netospheric origin  and  propagated  as  hydromagnetic  waves 

Axford  (1969)****  has  estimated  the  magnitude  of  three  forms  of  energy  dissi- 
pated during  a geomagnetic  substorm  to  be: 

1 ~1  o 18  erg/s  from  production  of  aurora 

2.  Joule  heating  (current)  of  polar  ionosphere  — 1 0 * ^ erg/s 

3.  Ring  current  10^  - 10^  erg/s 


*Akasofu,  SI,  “The  solar  wind  - magnetosphere  dynamo  and  the  magnetospheric  substorm,”  Planet 
Space  Sci,  23,  817, 1975 
**1  erg/s  = 0.1  „W 

***Akasofu,  SI,  “Polar  and  magnetospheric  substorms,”  Springer  Verlag,  1968 
****Axford,  I,  “Magnetospheric  convection,”  Rev  of  Geophys,  7,921,  1969 
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It  can  be  seen  that  the  total  rate  of  energy  release  is  on  the  order  of  the  energy  input  — 
indicating  that  substorms  must  be  frequent  phenomena,  assuming  substorms  are  the  major 
release  form. 

Each  substorm  injects  protons  into  the  trapping  or  radiation  belt.  The  adiabatic 
motion  of  these  protons  around  the  earth  forms  what  is  called  a "ring  current.”  If  there  has 
not  been  another  substorm  within  6 hours,  most  of  the  protons  will  have  been  lost  from  the 
ring  current.  However,  successive  substorms  at  intervals  of  2-3  hours  will  accumulate  pro- 
tons in  a very  intense  ring  current  (see  fig  1 ).  The  effect  of  this  ring  current  is  such  that  its 
magnetic  field  opposes  the  normal  geomagnetic  field,  and  so,  at  low  latitudes,  there  will  be 
a ground-level  geomagnetic  field  decrease  — called  a “geomagnetic  storm.”  The  develop- 
ment of  the  intense  ring  current  through  the  series  of  magnetospheric  substorms  occurring 
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Figure  1.  Relationship  between  substorm  occurrence  frequency  and  “trapped” 
proton  population  in  the  ring  current. 
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back  to  back  is  commonly  called  a magnetospheric  storm,  of  which  the  geomagnetic  storm  is 
one  manifestation 

The  rest  of  this  report  covers  in  detail  the  phenomena  of  magnetospheric  substorms. 
Succeeding  sections  provide  a refresher  on  magnetospheric  configuration  and  morphology; 
deal  with  the  phenomenology  of  the  "normal”  substorm;  give  some  experimental  observa- 
tions during  substorms;  and  examine  some  of  the  present  areas  of  conflict  in  the  interpreta- 
tions of  the  physics  of  substorms. 


MAGNETOSPHERIC  CONFIGURATION  AND  MORPHOLOGY 

The  solar  wind  is  a plasma  streaming  radially  outward  from  the  sun.  with  a weak  mag 
netic  field  imbedded.  Because  the  proton  gyro  radius  is  on  the  order  of  100  km  (near  the 
earth),  the  earth’s  dipolar  magnetic  field  acts  as  an  obstacle  to  the  solar  wind  flow  in  a basi- 
cally hydrodynamic  interaction.  The  following  descriptions  are  depicted  in  figures  2 and  3. 


MAGNETOPAUSE 


NOON  SOUTH 

boundary  or 
oval.  <3 


TROUGH 


magnetopause 

CURRENTS 


OUTER  LOOP  CURRENT 


MIDNIGHT  SOUTH 
BOUNDARY  or 
OVAL 


boundary  layer 


plasmasphere 


Figure  2.  The  magnetosphere 


Because  the  earth’s  dipolar  magnetic  pressure  falls  off  as  r . at  some  point  the  solar 
wind  pressure  (p)  dominates  the  earth's  magnetic  pressure  (ie.  p > B-/2p0)  - this  is  the  solar 
wind  region.  Closer  to  the  earth,  the  geomagnetic  pressure  dominates  (p  < B~/2p0)  and  we 
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f igure  3.  Magnetosphere,  showing  magnetic  configuration  and  bowshock. 

are  in  the  magnetosphere.  The  boundary  between  these  conditions  (where  p - B ~/2p0)  in 
called  the  magnetopause , and  moves  with  changes  in  the  dynamic  solar  wind  pressure. 

The  solar  wind  is  supersonic  (specifically,  super  Alfvenic)  and.  therefore,  a shock 
front  must  develop  in  front  of  the  magnetopause  (see  fig  3).  Between  the  shock  front  and 
the  magnetopause  is  a region  of  "shocked"  interplanetary  plasma  called  the  magnetosheath. 
In  this  region  are  "scrambled"  magnetic  field  lines,  which  must  finally  obey  the  magneto- 
pause boundary  condition  that  Bnorma|  disappear. 

The  supersonic  How  also  ensures  that  a cavity  of  extremely  low  plasma  density 
develops  b;  fi  nd  the  earth.  The  cavity  behind  the  earth,  caused  by  the  fluid  flow,  is  not  an 
open  cavity  because  of  the  isotropic  thermal  pressure  tending  to  close  it.  This  magneto- 
spheric  tail  is  complicated,  and  cannot  be  explained  by  simple  hydrodynamic  theory.  It 
must  be  treated  as  current  sheets  and  changing  magnetic  fields.  The  tail  field  lines  connect 
to  the  polar  regions  on  the  earth. 

A way  to  understand  the  dynamics  of  the  magnetosphere  is  to  use  the  concept  of 
merging  (reconnection)  of  B fields  imbedded  in  plasmas.  As  the  inner  field  lines  are  being 
pressed  together,  they  reconnect  to  form  a new  pair  (fig  4).  in  an  attempt  to  lower  the  total 
energy  in  the  field.  The  plasma  at  the  reconnection  point  (the  so-called  neutral  point)  is 
squeezed  and  accelerated  away  from  the  neutral  point.  The  exact  reconnection  process  is 
as  yet  poorly  understood. 

The  major  neutral  point  in  the  magnetosphere  occurs  where  the  interplanetary  field 
merges  with  the  dayside  outermost  closed  terrestrial  lines  (fig  5).  The  motion  of  the  solar 
wind  then  "drags”  these  lines  back  into  the  magnetotail.  A southward  interplanetary  mag- 
netic field  (IMF)  will  cause  enhanced  merging  at  the  front.  That  a northward  IMF'  can  have 
merging  is  one  of  the  subtleties  of  the  merging  theory,  because  the  fields  are  not  strictly  in 
antialignment.  The  effect  of  the  southward  IMF  is  that  more  field  lines  are  moved  across 
into  the  tail  in  a given  time,  causing  an  expansion  of  the  polar  cap/auroral  oval,  as  well  as 
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point. 

an  equatorward  motion  of  the  noon  elet't  (the  funnel-shaped  connection  between  the  niag- 
netosheath  and  the  earth,  tig  6). 


, Figure  6.  Fquatorward  motion  of  the  polar  cusp  (cleft)  due  to  southward 

turning  of  the  IMF. 

. The  increased  density  of  field  lines  in  the  tail  may  also  cause  merging  there  - far 

down  the  tail  (>  50  R|.  ).  This  tail  merging  will  be  a slow,  quasi-adiabatic  effect,  not  the 
rapid  substorm-caused  merging  to  be  examined  later. 

The  solar  wind  "blowing"  across  the  open  field  lines  (polar)  gives  a Vsoj  wmj  X Bn 
electric  field,  yielding  a 40-50-kV  potential  across  the  dawn-dusk  magnetotail.  This  action 
can  be  viewed  as  a dynamo,  with  a resulting  “crosstail”  current  (lowing  (fig  7).  Because  of 
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Figure  7.  Cross! ail  and  Held  aligned  currents  in  the 
magnetosphere. 

large  conductivities  in  the  nightsidc  auroral  oval,  there  may  be  slight  leakage  currents  (field 
aligned)  through  the  auroral  zone. 

SUBSTORMS 

I 

CLASSICAL  SUBSTORM 

After  a southward  turning  of  the  IMF.  there  is  an  increased  rate  of  merging  at  the 
magnetospheric  nose,  causing  the  neutral  point  to  move  earthward.  At  the  same  time,  the 
polar  cusp  moves  equatorward  and  the  auroral  zone  expands  (fig  (>).  The  magnetic  flux 
added  to  the  tail  from  the  solar  wind  increases  the  magnetic  field,  causing  an  increased  dawn/ 
dusk  crosstail  E-field  (there  is  also  an  increase  in  the  E-field  across  the  polar  cap).  The  in- 
creased tail  E-field  causes  an  H x B convection  of  the  tail  plasma  earthward,  which  creates  a 
thinning  of  the  tail  plasma  sheet.  This  enhances  the  reconnection  of  the  tail  x-type  neutral 
point  which,  in  turn,  causes  increased  (low  in  the  earthward  and  antiearthward  directions 
(fig  8). 

Svalgaard  ( I ‘>73)*  sees  this  as  a “local  collapse"  or  disruption  in  the  crosstail  current 
because  suddenly  there  are  no  particles  to  carry  the  current.  Hence,  the  disrupted  crosstail 
manifests  itself  in  a new  circuit:  from  the  dawn  side  magnetotail  down  the  field  to  the  dawn 
side  of  the  oval,  around  the  nightsidc  oval  (and  over  the  polar  cap),  finally  in  field-aligned 
How  up  to  the  dusk  side  tail  (fig  7). 


• Svalgaard,  L,  "Geomagnetic  responses  to  the  solar  wind  and  solar  activity,"  Stanford  University,  Institute 
for  Plasma  Research,  Rept  555,  1973 


Iiguic  S.  ('lunges  m the  magnetos  phene  conllgmalton  due  to 
southwaid  turning  of  the  IMI 

Akasolu  ( I ‘>70)*  has  a different  view  of  the  triggering  mechanism,  with  the  following 
time  sei|uenee 

1 I he  increased  crosstail  potential  causes  increased  current  flow  through  the 
auroral  circuit 

2 I he  increased  field-aligned  current  How  passes  some  threshold,  allowing  sudden 
enhancement  m anomalous  resistivity  along  a lie  let  line  (Akasolu  suggests  an  ion-acoustic 
wave  mechanism,  by  which  energy  is  taken  from  the  particle  How  and  converted  to  acoustic 
energy ) 

.)  1 lie  potential  drop  along  the  field  line  increases,  causing 

4 Acceleration  ol  auroral  particles  down  field  line  (f -field),  creating 

5 Increased  I region  ioni/ation  resulting  in  enhanced  conductivity  in  the  auroral 

/one. 

t>  I lie  enhanced  polar  circuit  conductivity  acts  as  a short  circuit  to  the  magneto- 
spheric  dynamo  and  disrupts  the  crosstail  current,  causing 

7.  f urther  increases  in  polar  circuit  current,  which  in  turn 

8 Increases  the  production  of  ion  acoustic  waves,  maintaining  potential  drop  along 
field  line  (maybe  even  increasing  it)  so  more  particles  are  accelerated,  creating  higliei  con- 
ductivity in  the  auroral  I region,  with  the  cycle  from  5-8  developing  in  a positive  feedback 
mechanism 


Akasolu  stales  that  the  short-circuiting  of  part  of  the  magnclotail  current  (through 
the  ionosphere)  w ill  cause  the  formation  of  a new  neutral  point  at  the  weakened  location  in 
the  magnclotail  Notice  that  the  current  flow  from  the  dawn  side  to  dusk  is  really  mani- 
fested by  electrons  flowing  the  other  direction,  ie.  streaming  down  the  nighlsulc  field  lines 


* Akasolu.  SI,  "Magnetospheue  substonus  a model,"  Solai  lmestiul  Physics/ 1 **  70,  pail  III,  p I HIM. 

Ox  lord  l*ie*».  l*>7  ’ 
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RECOVERY  PHASE 


Akasofu  yives  two  possible  causes  for  the  onset  of  the  recovery  phase: 

1 . A change  in  the  initial  external  stress  (a  turning  of  the  IME  from  southward  to 
northward). 

2.  The  ionosphere  may  be  able  to  hold  stress  only  temporarily  because  a neutral 
wind  is  generated  by  the  electrojet  in  the  auroral  E-region;  and  when  the  speed  becomes 
comparable  to  the  E-region  ion  speed,  the  electrojet  ceases  flowing  because  of  a V()  X B po- 
tential resisting  the  electrojet  motion. 

The  rest  of  the  recovery  phase  is  the  process  of  transporting  the  magnetic  flux  pro- 
duced during  the  storm  back  to  the  dayside  from  the  tail. 


RING  CURRENT  AND  TRAPPED  PARTICLES 

Many  of  the  particles  accelerated  earthward  by  the  tail  neutral  point  find  themselves 
on  "closed"  geomagnetic  field  lines  and,  once  there,  cannot  easily  diffuse  across  the  closed 
lines.  The  particles  in  the  "trapping  /.one”  undergo  drifts  around  the  earth  (fig  9)  caused  by 
the  curvature  and  gradients  in  the  earth’s  magnetic  field.  The  electrons  move  towards  the 
dawn  side  (eastward)  while  protons  drift  towards  evening  (westward),  constituting  a net 
westward  ring  current  in  the  trapping  zone.  The  trapped  electrons  have  the  same  energy 
(hence,  greater  velocity)  that  the  trapped  protons  have  (they  have  undergone  the  same 
accelerative  forces  in  the  tail  neutral  point)  and,  therefore,  display  a steeper  pitch  angle  dis- 
tribution; ie,  the  field-aligned  momentum  is  greater.  This  means  they  penetrate  deeper  into 
the  atmosphere  and  are  more  likely  to  undergo  collisions  that  scatter  them  out  of  the 
“trapped”  population.  Therefore,  mainly  protons  are  left  in  the  ring  current,  and  because 
these  protons  have  drifted  to  the  evening  sector,  there  remains  a net  asymmetry  in  the  ring 
current. 


figure  9.  Development  of  the  ring 
current. 


OBSERVATIONS  OF  SUBSTORM  ACTIVITY 


AURORAL  ACTIVITY 

Hie  close  relationship  between  auroral  displays  aiul  geomagnetic  disturbances  has 
been  recognized  since  ll>5l).  Akasofu  ( l%8)  compared  aurora  to  the  cathode-ray  tube;  ie, 
the  tluorescent  display  as  the  screen  is  bombarded  by  the  electron  beam.  The  aurora  are 
large-scale  phenomena  with  complicated  spatial  features  that  vary  in  time,  Each  observing 
station  travels  through  the  aurora  as  the  earth  rotates,  making  it  difficult  to  understand  the 
auroral  morphology  without  the  use  of  all-sky  camera  chains  and  polar  satellites.  Akasofu 
t P>o4)*  has  developed  a schematic  chain  of  events: 

I T*0  quiet  (fig  1 0 A ) 

2.  1=0-5  min  brightening  of  the  quiet  midnight  arc  (fig  1 OB) 

.V  1=5-10  min  poleward  motion  of  the  brightened  arc,  trailing  “breakup” 

(bulge  effect)  (fig  IOC  ) 

4 1 = 10-30  min  bulge  expands  mainly  through  a westward  traveling  surge 

(fig  1 01)) 

5.  T*30  min- 1 h recovery  stage:  contraction  of  the  poleward  bulge  (fig  1 OF) 

The  initial  poleward  motion  of  the  aurora  (step  3)  can  be  explained  by  the  action  of 
enhanced  reconnection  in  the  tail  (fig  I I ).  because  auroras  tend  to  appear  on  the  poleward 
edge  of  the  close  field  lines.  The  westward  traveling  surge  seen  in  the  aurora  (step  4)  is 
caused  by  an  influx  of  energetic  electrons  the  auroral  motion  is  a manifestation  of  the 
westward  motion  of  the  particle  dumping  region  Later  in  this  report,  in  examining 
Akasofu’s  models  (fig  I 3),  the  substorm  cause  of  the  westward  motion  will  be  pointed  out. 

SATELLITE  MEASUREMENTS 

Satellite  measurements  have  given  a picture  of  the  magnetospheric  configuration 
during  substorms.  There  are  basic  problems  in  interpreting  the  information,  however,  be- 
cause it  is  still  difficult  to  separate  space  and  time  variations. 

During  the  storm  onset  there  is  plasma  motion  in  the  tail,  indicating  a thinning  of  the 
plasma  sheet  and  earthward  motion  of  the  trapping  regions  (see  fig  8).  This  occurs  because 
the  increased  convection  causes  a tail  E-field  pinch.  The  geomagnetic  field  configuration  at 
6-10  K|.  becomes  more  dipole-like  coinciding  with  westward  electrojet  m the  auroral 
region.  This  increased  activity  causes  the  injection  of  plasma  into  the  trapping  regions  from 
the  tail  regions. 

SUBSTORM-RELATED  CURRENT  SYSTEMS 

Chapman  and  Birkeland  in  I ‘MO  developed  “equivalent  current"  systems  to  describe 
ground  magnetic  fluctuations,  in  which  they  observed  the  need  for  an  auroral  electrojet 

•Akasofu,  SI,  “The  development  of  the  auroral  substorm,"  Plant  Spaa1  Set,  l_2,  273,  I'tp-t 
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Figure  10.  Development  of  the  auroral  substorm  (Akasofu,  1 064,  reprinted  with  permission). 
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Figure  1 1 . Poleward  motion  of  (he  aurora  due  to 
enhanced  reconnection  in  the  niagnetotail. 


circuit,  along  with  the  field-aligned  (or  Birkeland)  currents.  The  following  models  of 
Akasofu  (and  others)  use  current  systems  consisting  of  day  /evening  asymmetric  ring  current, 
field-aligned  currents  into  the  dawn/dusk  oval,  and  the  auroral  electrojet  (fig  I 2).  The  low- 
latitude  magnetic  field  measurements  require  an  eastern  (dusk/dawn)  tail  current  (opposing 
the  normal  dawn/dusk  current).  Akasofu  recognized  that  the  loss  of  westward  current  is 
equivalent  to  the  needed  eastward  current,  and  so  he  allowed  the  tail  current  to  he  shunted 
into  polar  current  (fig  1 2 A).  Since  this  recognition,  people  have  worked  on  more  sophis- 
ticated 3D  models  (fig  I 2B)  (Fukushima.  1972;  Yasuhara  and  others)*  that  include  the 
partial  ring  current 

Akasofu.  in  his  early  work  ( 1970).  also  developed  a concept  of  the  simple  time  vari- 
ations of  the  current  system  during  a substorm.  Notice  from  figure  1 3 that  the  region  of 
reduced  plasma  density  spreads  radially  outward  along  the  neutral  sheet.  The  spreading 
down  the  tail  cause's  the  poleward  motion  of  the  aurora,  explained  earlier,  while  the  cross- 
tail  spreading  causes  the  endpoints  of  the  field-aligned  magnetospheric  currents  to  move 
westward,  creating  a westward  surge  of  auroral  activity. 

AREAS  OF  CONFLICT:  AKASOFU  AND  THE  DYNAMIC  AURORAL  OVAL 

Akasofu  first  pointed  out  in  1973  that  the  oval’s  position  varies  with  the  direction 
of  the  interplanetary  field,  thus  affecting  auroral  electrojet  ( Al ■')  measurements. 


*i  ukushima.  N,  "Polai  magnetic  suhstomis."  Planet  Space  So.  20.  1443,  l‘>72 
N asuhara.  f . K amide.  Y.  and  Akasofu.  SI.  "A  model  of  the  magnetospheric  substorm,"  Planet  Space  Sci. 
23,  S4S,  |<>75 
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Figure  I2A.  Substomi  current  systems  after  Akasofu  ( ll>72 ). 


Figure  13.  Development  of  a niagnetosphcrie  substorm, 
as  the  crosstail  current  is  deflected  to  the  polar  region. 
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In  a process  exactly  like  that  indicated  in  the  early  development  of  substorms 
(McPherron,  1072*).  the  enhanced  merging  at  the  magnetospheric  nose  due  to  the  south- 
ward turning  of  the  IMF  moves  the  polar  cusp  cquatorward  and  causes  an  expansion  of 
the  open  polar  cap  (auroral  oval)  (see  fig  6,  1 1 ). 

Akasofu  contends  that  variations  in  this  near-reversible  process  are  continually 
present,  and  are  the  cause  of  many  of  the  conflicts  that  will  be  mentioned  later.  Specifi- 
cally. motions  of  the  oval  (and  its  always-present  but  slight  electrojet)  will  cause  ground- 
level  magnetic  variations  that  may  mask  or  obscure  the  actual  substorm  onset  as  observed 
by  ground-level  magnetometers.  More  important,  however,  is  the  fact  that  when  the  IMF 
has  a large  northward  component,  the  nose  merging  decreases  to  near  zero,  and  the  auroral 
oval  contracts  far  poleward  (up  to  ~80°).  If  the  oval  is  above  ~70°  dipole  latitude,  the 
magnetic  activities  associated  with  substorms  will  not  be  seen. 

Therefore,  as  Akasofu  has  been  strongly  pointing  out,  studies  that  use  the  auroral 
electrojet  indices  (or  magnetic  signatures  at  nominal  auroral  stations)  will  not  see  activity 
when  IMF  B.  is  northward,  and  so  all  "activity”  will  be  strongly  correlated  to  southward 
IMF. 


THE  GROWTH  STAGE:  AKASOFU  VS  McPHERRON 

The  basic  phenomenology  of  the  magnetospheric  changes  following  a southward 
turning  of  the  IMF  (given  by  McPherron  ( I ‘>72),  fig  8)  is  accepted  by  Akasofu  It  is 
MePherron's  interpretation  of  succeeding  events  that  Akasofu  takes  exception  to. 

McPherron  envisions  that  towards  the  end  of  the  "growth  stage,”  the  plasma  sheet  at 
~ i 0 R |.  becomes  thinner  (narrower)  because  of  increasing  tail  magnetic  pressures  caused  by 
a southward  IMF.  He  suggests  that  this  thinning  creates  a neutral  point  field  reconnection, 
marking  the  onset  of  the  expansive  stage  of  the  substorm.  In  this  theory,  the  tail  is  the 
location  of  the  triggering  instability. 

Akasofu,  on  the  other  hand,  believes  that  an  anomalous  resistivity  in  the  field- 
aligned  current  How  causes  enhanced  current  flows  along  the  field  lines  into  the  auroral 
oval.  This  earthward  llow  disrupts  the  crosstail  currents  because  there  is  no  plasma  to  carry 
it.  This  disruption  happens  at  the  same  tail  location  at  which  MePherron’s  model  displayed 
the  initial  thinning  (ie.  IORp).  The  tail  plasma  sheet  thins,  creating  a reconnection  point, 
and  the  substorm  expansion  phase  begins. 

The  argument  goes  deeper,  however.  As  Akasofu  stresses,  there  is  ample  evidence 
that  substorms  are  as  common  with  the  IMF'  northward  as  they  are  with  it  southward,  albeit 
the  substorm  energy  expenditures  are  way  down.  Remember  that  MePherron's  model  used 
the  enhanced  merging  from  the  IMF  being  directed  southward  to  trigger  the  substorm 
mechanism.  Unfortunately,  McPherron  deals  with  Akasofu 's  “northward"  events  by  ignoring 
them  on  the  grounds  that  the  total  energy  release  therein  is  small  and,  therefore,  “they  don’t 
represent  any  great  change  within  the  magnetosphere."  Akasofu 's  point  was  missed  entirely 
how  do  these  “northward”  substorms  trigger?  Certainly  enhanced  tail  pressures  cannot  be 


•McPherron,  Rl.,  “Substorm  related  changes  in  the  geomagnetic  tail: 
20,  1521,  I ‘>72 
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the  source,  and  then  McPherron  must  postulate  that  the  substorms  with  the  IMF  directed 
northward  are  intrinsically  different  from  the  substorms  with  the  IMF  southward,  contrary 
to  any  evidence  available  today. 

CONCLUSIONS 

The  basic  mechanism  of  magnetospheric  energy  release  through  substorms  is  be- 
coming understood  with  the  aid  of  satellite  measurements.  The  energy  dissipation  is  in  the 
auroral  E-region,  and  comes  about  through  particle  precipitation  and  current  (joule)  heating. 
This  energy  input  into  the  polar  ionosphere  and  neutral  atmosphere  is  the  cause  of  the 
ionospheric  substorm  portion  ot  Akasofu’s  substorm  manifestations,  and  the  ensuing  global 
ionospheric  disruptions  impact  upon  communication  systems  using  ionospheric  propagation. 
Therefore,  understanding  the  temporal  and  spatial  substorm  variations  should  aid  in  deter- 
mining the  temporal  and  spatial  variations  in  the  global  ionosphere.  Extending  the  present 
DMSP  satellite  capabilities  to  monitor  real-time  auroral  activity  would  provide  energy  con- 
centration inputs  for  ionospheric  storm  models.  Ionospheric  storm  modeling,  which  is  in 
its  early  stages  of  development,  is  discussed  in  chapter  2. 
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CHAPTER  2:  IONOSPHERIC  STORMS 


INTRODUCTION 

Ionospheric  storms  arc  the  F^-region  variations  (including  enhancements  and  deple- 
tions in  electron  density,  as  well  as  raising  and  lowering  of  the  Fs  peak  height)  associated 
with  magnetic  storm  variations  The  consensus  in  the  field  is  that  these  storm-time  iono- 
spheric changes  are  of  secondary  nature  caused  hy  the  geomagnetic  storm  aftereffects.  The 
resultant  changes  are  an  outcome  of  several  competing  and  cooperating  processes  and, 
the  a*  fore,  each  individual  stonn  pattern  appears  complex  and  unpredictable.  Several 
workers  have  developed  “average”  stonn  patterns  for  given  locations  as  well  as  general 
global-type  average  morphologies.  Unfortunately,  the  individual  stonns  do  not.  in  general, 
follow  these  average  events,  making  the  "average”  stonn  effectively  useless  for  practical 
propagation  predictions. 

The  major  problems  in  studying  ionospheric  stonns  stem  from  the  total  number  of 
variables  that  can  affect  the  electron  density  changes.  These  include  ( 1 ) the  time  of  the 
sudden  commencement  if  the  storm  starts  during  the  high-electron-density  day,  the 
effects  could  be  different  than  at  night;  (2)  the  state  of  the  neutral  thennosphere  this 
controls  the  whole  chemical  loss  process  and  may  affect  the  propagation  of  the  disturbance 
equatorward;  (3)  the  actual  magnetic  field  and  electrical  current  variations  these  are  the 
driving  force  behind  the  ionospheric  disturbance;  (4)  the  position  (on  the  earth)  of  the 
recording  station  relative  to  the  magnetic  and  electric  current  changes  a short-time  scale 
variation  will  cause  very  localized  effects  (Park.  1474).*  It  is  likely  the  neglect  of  one  (or 
more)  of  the  above  factors  in  every  analysis  has  led  to  the  apparent  inability  to  understand 
ionospheric  stonn  behavior. 

Ionospheric  stonns  are  observed  with  two  basic  experimental  probes:  ( 1 ) hf  vertical 
incidence  sounders  yielding  basically  the  density  and  altitude  of  the  F-region  peak;  and 
(2)  total  electron  content  (TEC)  systems  which  use  the  Faraday  rotation  of  vhf  satellite 
signals  to  give  the  integrated  electron  density  along  the  signal  path.  The  measurements  are 
complementary  because  the  TEC  system  is  also  affected  by  the  electron  density  above 
the  F-region  peak,  a region  which  the  vertical  incidence  sounder  can  obviously  not  probe. 

This  report  reviews  a large  set  of  the  F-region  ionospheric  storm  analyses,  the  iono- 
spheric processes  believed  to  be  involved,  and  the  "state  of  the  art”  in  storm  prediction. 

A final  section  investigates  the  usefulness  of  this  state  of  the  art  to  real-time  propagation 
prediction  and  hf  frequency  management.  It  also  suggests  future  study  to  enable  the 
definition  of  storm  morphology. 


•Park,  CG,  "A  morphological  study  of  substorm-associated  disturbances  in  the  ionosphere,”  J Geophys 
Rev.  79.  2821,  1974 
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FEATURES  OF  THE  IONOSPHERIC  STORM 


AVERAGE  STORM  EFFECTS 

Because  of  the  large  storm-to-storm  variations,  many  workers  have  developed 
“average  storm  patterns”  in  attempts  to  describe  and  understand  ionospheric  storm  behavior. 
Mendillo  ( 1971)*  (using  TEC  measurements)  observed  that  there  are  two  distinct  time 
patterns  superposed  on  these  storms:  ( 1 ) storm-time  variations  ( Dcy ) reckoned  from  time 
of  sudden  commencement;  (2)  local  time  variations  (SI))  reckoned  by  local  time  with  the 
sudden  commencement  in  day  1.  The  (l)jj  j-)  (storm  time)  variations  have  a distinct  pattern 
(fig  14A),  with  a 12-hour  positive  phase  followed  by  a 72-hour  negative  phase.  The  (SD) 
(local  time)  pattern  (fig  14B)  indicates  that  this  large  positive  effect  is  at  the  sunset  period 
on  the  first  day  of  the  storm.  Again,  notice  that  any  given  storm  may  not  follow  this 
average  behavior.  In  particular,  the  positive  phase  appears  to  be  dependent  on  the  time  of 
the  sudden  commencement,  and  many  storms  (usually  with  local  nighttime  sudden  com- 
mencements) have  no  positive  phase.  On  the  other  hand,  work  done  at  NOSC  indicates  that 
during  the  winter  season  the  whole  storm-time  effect  may  be  a positive  phase. 


HOURS  AFTER  STORM  SOI  SD2  SD3  SD4 

COMMENCEMENT  (SCI 

Figure  14.  TEC  “average  storm"  variations  showing  both  storm  time  ( l and  local  time  (SI)) 
dependence  (Mendillo,  M,  "Ionospheric  total  electron  content  behavior  during  geomagnetic  storms. 
Nature.  234,  23,  1971) 


Mendillo  states  that  this  analysis  was  used  to  look  for  seasonal  variations  in  storm 
behavior,  and  that  no  great  seasonal  differences  were  discovered  This  is  somewhat  surpris- 
ing because  the  above-mentioned  analysis  at  NOSC'  indicates  very  strong  seasonal  differences 
in  storm  variations. 


•Mendillo.  M,  "Ionospheric  total  electron  content  behavior  during  geomagnetic  storms."  Nature.  234. 
23, 1971 
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NEGATIVE  PHASE 


The  "average  storm-time  variation”  is  found  to  be  an  increase  in  electron  density 
(positive  phase)  followed  by  a decrease  (negative  phase).  In  this  section  the  negative  (or 
main)  phase  and  the  physical  mechanisms  are  discussed 

The  negative  phase  of  the  ionospheric  storm  appears  to  be  caused  by  the  neutral 
upper  atmosphere’s  reaction  to  the  energy  deposition  of  geomagnetic  storms.  This  energy 
is  deposited  in  the  upper  atmosphere  in  two  forms  particle  dumping  and  the  resistive 
"joule”  heating  of  electrical  currents.  The  particles  are  mainly  dumped  near  the  geomagnetic 
"cusp”  region,  although  there  is  a spread  around  the  auroral  zone  as  well  as  down  to  mid- 
latitudes. Recently  it  has  been  determined  that  the  electrical  currents  in  the  auroral  zone 
can  be  very  large  and  hence  cause  significant  joule  heating.  These  energy  inputs  modify 
the  high-altitude  global  wind  systems  and  also  cause  gravity  wave  trains  to  propagate  out- 
ward (equatorward  and  poleward)  with  the  apparent  effect  being  a turbulent  mixing  of  the 
neutral  atmosphere.  This  mixing  is  followed  by  an  increase  in  the  chemical  loss  rates  which 
should  result  in  a net  decrease  in  electron  density  as  well  as  an  altitude  increase  in  the  peak 
height. 

Other  explanations  for  the  electron  density  decrease  invoke  neutral  wind  or  electric- 
field  transport  mechanisms.  The  transport  to  the  higher-loss  region  necessitates  a downward 
motion  of  the  F-layer,  and  this  is  contradictory  to  experimental  evidence.  Therefore,  it  is 
usually  assumed  that  the  loss  rates  themselves  are  changed  during  the  negative  phase  of  the 
storm,  with  the  change  mechanism  being  the  neutral  thermosphere.  This  manner  of  change 
explains  both  the  decrease  in  electron  density  and  the  increase  in  the  altitude  of  the  F-region 
peak. 


POSITIVE  PHASE 

There  is  general  agreement  that  the  positive  early  storm  effects  are  caused  by  trans- 
port of  the  positive  ions  and  electrons.  Vertical  transport  into  lower-chemical-loss-rate 
regions  (ie,  higher  altitudes)  by  either  thermospheric  winds  or  electrodynamic  forces  is 
likely  to  dominate  any  horizontal  transport  mechanism  because  of  scale  height  effects. 

Prolss  and  Najita  ( 1975)*  suggest  that  the  electrodynamic  uplifting  is  responsible 
because  evidence  (both  experimental  and  theoretical)  is  that  the  wind  motions  would  be 
zonal  (east-west)  component  dominated  (causing  no  upward  transport  along  field  lines). 
Lanzerotti  and  others  (1975)**  showed  that  the  large  local  evening  enhancements  in  some 
storms  are  associated  with  positive  bay  geomagnetic  disturbances,  which  indicate  the  exist- 
ence of  large  electric  currents  that  would  create  upward  electrodynamic  forces.  Mendillo 
(private  communication)  has  observed  this  association  also,  but  finds  no  relationship  of  the 
electron  density  enhancement  magnitude  to  that  of  the  positive  bay.  Papiagannis.  Mendillo. 
and  Klobuchar  ( 1971  )***  found  that  the  maximum  total  electron  content  (TFT)  during 


•Prolss,  GW,  and  Najita.  K,  "Magnetic  storm  associated  changes  in  electron  content  at  low  latitudes," 

J Atmos  Terr  Phys,  37,635, 1975 

••Lanzerotti,  LJ,  Cogger,  LL,  and  Mendillo,  M,  "Latitude  dependence  of  ionospheric  total  electron 
content:  Observations  during  sudden  commencement  storms,”  J Gcophys  Res,  80,  1287,  1975 
•••Papagiannis,  MD,  Mendillo,  M,  and  Klobuchar.  JA,  “Simultaneous  storm-time  increases  of  (he  ionospheric- 
total  electron  content  and  the  geomagnetic  field  in  the  dusk  sector,”  Planet  Space  Sci,  ]9,  505,  1971 
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these  positive  phase  periods  averaged  much  higher  than  usual  local  noontime  value,  with  the 
positive  phase  being  a sunset  period  phenomenon 

LARGE  AND  R APID  DROPS  IN  TEC  NEAR  SUNSET  PERIOD 

In  examining  the  lit  data  of  AECKL’s  Sagamore  Mill  Radio  Observatory  at  Hamilton, 
Mendiilo  and  others  1 1^74  )*  observed  that  on  a few  occasions  an  abnormally  large  and  rapid 
drop  in  1 1 1 is  observed  near  the  sunset  period  They  find  the  decay  rates  tdN  | dt)  to  be 
more  than  IU  times  greater  than  the  respective  monthly  mean  Because  they  did  not  believe 
that  any  change  in  the  loss  process  could  be  k such  drastic  nature,  their  conclusion  was  that 
due  to  plasmasphenc  contraction  the  radio  propagation  path  was  crossed  by  the  midlatitude 
trough  Vlendillo  and  others  developed  a strong  case  lor  the  equatorward  movement  of  the 
plasmapause  boundary  causing  the  sharp  decreases  The  Hamilton  TI  C'  subionospheric 
point  l.  value  is  2.8,  and  movements  of  the  boundary  into  this  region  are  not  impossible 
However,  data  from  the  8 March  107Q  event  indicate  that  the  sudden  decreases  exist  down 
to  Arecibo,  with  an  1 value  of  I 4 (fig  15)  This  is  far  equatorward  of  any  reasonable 
plasmapause  boundary . and  so  the  Mendiilo  explanation  would  seem  unreasonable. 

Davies  ( ll>74)**  developed  an  explanation  for  the  sudden  decrease,  in  which  he 
allows  a sudden  atmospheric  mixing  to  cataclysmically  increase  the  loss  rates.  The  sudden 
mixing  is  caused  by  a “shock  front"  (see  appendix)  of  turbulence  propagating  equatorwards 
from  a "supersonic"  disturbance  in  the  auroral  /one.  I or  simplicity . Davies  considers  the 
disturbance  to  be  particles  dumped  out  of  the  magnetosphere  into  the  geomagnetic  cleft 
(cusp),  which  is  at  approximately  the  local  noon  longitude,  l or  propagation  velocities 
of  4°-5°  h ('•150  m s),  the  front  will  propagate  equatorwards  in  such  a manner  as  to  cause 
the  sudden  depletions  with  the  time  lags  derived  for  the  8 March  D>70  event  (see  appendix) 

REVIEW  OF  PAST  WORK  IN  IONOSPHERIC  STORMS 

SOME  OBSERVATIONS  AND  ATTEMPTS  AT  GLOBAL  MORPHOLOGY 

Many  workers  in  the  past  have  examined  data  from  one  or  more  stations  and 
developed  storm  morphologies.  One  of  the  first  was  a study  done  by  Matsushita  (1050)***, 
who  found  that  the  equator  region  had  positive  storm  effects  while  mid  and  high  latitudes 


•Men  till*'.  M,  and  klobuchar,  J,  An  Atlas  of  the  midlalitudc  F-region  response  lo  geomagnetic  storms. 
Al  l R1  Technical  Report  74-OOf>5,  February  1074 
••Davies,  k.  "Studies  of  ionospheric  storms  using  a simple  model,"  J Geopliys  Res,  TO.  t>05,  1074 
•••Matsushita.  S,  "A  study  of  the  morphology  of  ionospheric  storms,"  J (leoplivs  Res,  r>4.  505.  1050 
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COLUMNAR  ELECTRON  CONTENT,  IN  UNITS  OF  1017  ELECTRONS  m' 
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Figure  15,  TEC  variations  for  ionospiicric  storm  of  8 March  1*770  for  several 
stations.  The  dotted  lines  indicate  the  average  diurnal  patterns  for  each  station. 
(Klobuchar,  JA,  and  others,  “Ionospheric  storm  of  March  8,  1*770,"  J Gcophys  Res, 
76,6202.  1*771) 
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luul  negative  effects  More  recent  works  (Prolss  aiul  Nuiita.  1 075 . Kane.  1^7.*.  Kane. 

I‘>75;  Kush,  1 7 2 > * have  “typical''  storm  patterns,  hut  to  quote  Kane  ( 1073).  "deviations 
from  the  average  pattern  seem  to  he  mom  a rule  than  an  exception."  I he  basic  pattern 
appears  to  he  a positive  phase  on  clay  I . followed  hy  a negative  storm  phase  on  day  2. 
recovering  somewhat  on  day  3 Metulillo  ( I ‘*7 1)  and  Jones  ( l‘>7 1 )•*  indicate  that  if  the 
sudden  commencement  occurs  before  sunrise,  them  will  probably  be  no  positive  phase 
Kane  ( N73)  found  that  at  low  latitudes  the  “average  pattern"  of  lt,l  > increases  is  rarely 
seen  often  them  are  large  decreases 

Both  Kane  ( 1073)  and  Kush  ( I ‘>72)  explain  that  often  different  ionospheric  behavior 
patterns  are  observed  between  locations  separated  by  a few  thousand  km.  Kane  postulates 
that  ionospheric  storms  are  the  result  of  local  "typhoon"  type  events  in  the  thermosphere, 
and  hence  we  need  a flo.se  global  network  ol  stations  monitoring  f(,T  > and  l i t'  These 
disturbances  "originate  in  polar  and  auroral  regions  and  forming  convection  cells,  descend 
to  lower  latitudes  in  turbulent  fashion  " Hanks  and  Nagy  ( I ‘>741***  also  forward  this 
cyclonic  theory  of  ionospheric  disruptions  m order  to  explain  the  localized  effects  some 
ionospheric  storms  cause 

Park  ( |074)  described  an  observation  of  an  isolated  geomagnetic  substorm  of 
moderate  intensity  tKj,  401  that  produced  a /< »« ionospheric  storm  over  the  eastern 
continental  United  States  In  that  local  area  this  ionospheric  disturbance  was  similar  to  the 
behavior  seen  during  major  ionospheric  storms  1 1 is  conclusion  is  that  "ionospheric  storms 
can  be  understood  in  terms  of  the  superposed  effects  of  substorms."  The  impact  of  tins  is 
that  the  planetary  geomagnetic  disturbance  indices  would  not  be  expected  to  relate  strongly 
to  many  ionospheric  disturbances  in  a given  localized  area,  and  that  mom  “zonal"  parameters 
must  be  used  to  understand  ionospheric  behavior 

TIUOKI  riCAl.  TUI  KMOSPHTKIC  K1SPONSI 

Kiclunotul  and  Matsushita  t 1075)****  developed  a computet  model  of  the  thermo 
spheric  response  to  magnetic  substorms  They  conclude  that  the  generation  of  large  scale 
gtavity  waves  causes  an  impulse-like  disturbance  front  to  propagate  poleward  and  equator 
waul  at  750  m/s  The  amplitude  of  the  front  decays  very  little  before  reaching  the  equator 
Also,  ion  drag  winds  am  set  into  motion  at  the  auroral  oval,  but  because  no  diurnal  variations 
am  included  in  the  model,  them  am  no  longitudinal  effects. 

— 

•Kane.  KP.  "Storm  time  variation  of  T.'  region."  Annls  (leophvs,  2‘>.  2'.  I *•> 7 A 
Kane.  KP.  "(ilohal  evolution  ol  I he  ionospheric  electron  content  dining  some  magnetic  slot  ms." 

J Atmos  and  Terr  Plivs.  37.  nOI . I *>75 

Kush.  CM.  "Some  effects  of  neutral  wind  changes  on  the  low  latitude  T regton,"  .1  Atmos  Ten  I’hvs, 

34.  I40t.  I»72 

••Jones,  Kl  . "Storm  tune  variation  ol  TJMnyet  electron  concentration,"  J Atmos  and  Ten  I’hvs,  tj. 

37*).  I *171 

•••Hanks,  PM.  and  Nagy.  AT.  "Cyclonic  disturbances  and  llreii  consequences  In  the  thermosphere," 
(Jeophys  Kes  letters,  1.305,  l‘>75 

••••Klchmond.  Al),  and  Matsushita.  S.  "Thermospheric  response  to  a magnetic  suhsloim,"  J (ieophvs 
Kes.  HO.  2HJ0,  P>75 
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Notice  that  the  750-m/s  disturbance  propagation  velocity  is  considerably  higher  than 
the  150  m/s  inferred  from  data.  It  may  be  that  the  calculated  winds  following  the  disturbance 
with  velocities  of  200  m/s  are  the  mixing  source. 


RESULTS  OF  PRESENT  WORK 

We  have  taken  vertical  incidence  sounder  data  from  Vandenberg,  California,  and 
Honolulu,  Hawaii,  for  several  ionospheric  disturbances  between  May  1973  and  the  present. 
Using  a pre-event  average  defined  by  the  5 days  prior  to  the  sudden  commencement,  we 
calculated  the  average  variations  of  the  storms  as  Mendillo  did.  We  found  a distinct  seasonal 
difference  that  was  separated  as  summer  (May-Aug)  and  winter  (Nov-Feb)  (fig  16A,  B). 

A.  HAWAII  SUMMER 


HAWAII  WINTER 


Figure  16A,  B.  Average  fQF2  storm  variations  (using  local  time)  for 
Hawaii,  summer  and  winter  storms. 

The  Hawaii  station  quiet  average  day  had  a very  peaky  density  function  (around  noon) 
for  most  events.  This  fact  made  the  Honolulu  data  rather  intractable  because  slight  shifts 
(along  the  time  axis)  of  the  storm  data  would  cause  large  enhancements  and  depletions  even 
if  the  daily  maximums  did  not  vary.  These  variations,  when  averaged,  almost  cancel  out  and 
yet  leave  large  standard  deviations  in  the  error  analysis.  Basically,  the  problem  is  that  the 
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ionospheric  stonn  modified  the  diurnal  variations  in  a presently  unpredictable  manner  In 
an  attempt  to  find  some  predictable  facet  of  a storm,  we  next  took  the  daily  peak  (ignoring 
time  information)  and  developed  an  average  storm  behavior  for  Hamilton  (Wallops  Island 
AFCRL  data),  Vandenberg,  and  Hawaii  (fig  1 7A.  B,  C).  Again,  a strong  seasonal  dependence 
is  noted.  Notice  that  Wallops  Island  summer  shows  a large  increase  on  day  1 followed  by 
a strong  decrease  on  day  2,  with  a recovery  phase  on  days  3 and  4.  Ixannning  the  second 
day  of  the  stonn,  one  discovers  that  four  of  the  24  summer  events  studied  in  Wallops 
Island  data  actually  show  increases  above  nonnal  “quiet  day  ." 
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Figure  17.  Storm  FqF-,  variations  (daily  maximum)  for  summer,  winter,  and 
equinoctial  storms. 

As  another  exercise,  the  Wallops  Island  data  were  used  to  compare  day  3 to  day  2. 

In  14  out  of  26  summer  events  the  fQF  ->  value  on  day  3 was  higher  than  on  day  2;  so  far 
we  have  been  unable  to  separate  those  14  from  the  remaining  I 2 by  using  any  other  para- 
meter. This  means  that  on  day  2 of  the  event,  present  information  does  not  allow  prediction 
on  whether  day  3 will  be  more  or  less  severe  than  day  2. 

The  above  indicates  that  although  an  "average"  stonn  may  be  defined,  we  cannot 
accurately  estimate  the  ionosphere  variations  for  a given  event  by  using  the  presently 
available  information.  Although  a storm  generally  shows  decreases  in  the  F-region  electron 
density,  neither  the  amount  of  this  decrease  nor  its  temporal  variation  appears  related  to  the 
magnetic  stonn  variations  as  described  by  Ap  values. 
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EVALUATION  OF  AN  OPERATIONAL  MODEL 


Patton  and  Cohen  ( 1971  )*  have  developed  a prediction  model  in  order  to  allow  the 
USAF  to  use  real-time  frequency  management  on  certain  hf  systems  (OTH  radar).  These 
authors  decided  to  deal  with  only  negative  storm  phases  (and  predictions),  which  have  the 

* largest  system  impact  because  they  cause  hf  system  outages.  The  model  was  developed  to 
fit  the  storm  morphology  described  in  Obayashi  (1964)  and  Davies  (1965)  (Davies  tabulated 
Matsushita’s  ( 1959)  data  on  average  behavior).**  Davies  states  of  his  table:  “the  overall 
average  is  therefore  not  representative  of  a single  storm  in  either  of  those  seasons  (summer, 
winter).”  The  data  indicate  an  approximate  15%  decrease  at  the  equator  varying  to  ~40% 
decrease  at  60°  latitude. 

The  model  incorporates  seasonal  and  latitudinal  variations  by  means  of  a “local  K 
index”  Kl  which  is  averaged  with  the  Kp:  K = ( 1 .5Kp  + Kl)/2.5.  Kl  has  a maximum  at 
summer  high  latitude  and  a minimum  at  winter  low  latitude  in  the  form: 
summer  winter 

0°  KL  = 3 Kl=2 

90°  KL=9  Kl=7 

The  storm  MUF  is  then  described  by: 

MUF  storm  _ 13.1  - K 

MUF  quiet  11.1 

to 

This  model  predicts  average  ionospheric  storm  responses  according  to  Obayashi  and 
Matsushita  morphologies.  As  described  in  the  last  section,  the  data  studied  at  NOSC  indicate 

* that  lower-latitude  regions  (Hawaii)  may  be  considerably  more  responsive  than  midlatitude 
regions.  Also,  in  18  winter  events  measured  at  Wallops  Island,  Virginia,  only  one  decreased 
to  significantly  below  average.  In  other  words,  the  above-described  model  does  not  match 
the  ionospheric  storm  responses  accumulated  at  NOSC. 

FUTURE  CONSIDERATIONS 

It  has  been  emphasized  that  the  presently  available  data  set  is  not  spatially  dense 
enough  to  develop  an  understanding  of  ionospheric  storms.  The  future  should  thus  initially 
include  the  operation  of  such  a measurement  network  covering  at  least  the  Pacific  area  so 
that  complete  storm  morphologies  can  be  developed.  Once  this  assessment  process  is 
completed,  then  perhaps  a prediction  capability  will  be  at  hand.  A stop  gap  measure  being 
undertaken  by  Blei weiss  and  Argo  (private  communication)  is  the  using  of  the  AWS  sounder 
network  surrounding  the  Pacific  in  a real-time  ionospheric  assessment  program.  We  now 

* know  that  the  wide  separation  of  the  sounders  will  probably  allow  only  very  degraded 
analysis,  and  so  further  effort  must  be  expended. 


•Patton,  DE,  and  Cohen,  ML,  “A  simple  magnetic  storm  prediction  model  suitable  for  real-time  frequency 
management,”  Raytheon  report  published  at  ARP AOHD  Tech  Review,  17-18  March  1971 
**Obayashi,  T,  Research  in  Geophys,  chapter  14,  vol  1 , MIT  Press,  Cambridge,  1964 
Davies,  K,  Ionospheric  Radio  Waves,  Blaisdale  Publishing  Co,  Walthrow,  MA,  1965 
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In  short,  the  ionospheric  storm  can  be  viewed  as  a meteorological  process  (albeit 
thermospheric)  with  the  attendant  prediction  problems.  A close  network  of  sounders  will 
allow  the  measurement  of  “frontal"  motion  and  the  movement  of  “high”  and  “low”  density 
regions.  But  until  these  highs  and  lows  can  be  viewed  as  finite,  moving  regions,  their 
effects  on  a single  ionospheric  measuring  system  must  appear  capricious. 
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APPENDIX: 

GEOMAGNETIC  STORM  EFFECTS  ON  THE  MIDLATITUDE  IONOSPHERE: 

A "SHOCK  FRONT”  TURBULENT  MIXING  MODEL  AND 
CALCULATIONS  OF  THE  FRONT  SHAPE 

INTRODUCTION 

Ionospheric  storms  are  commonly  caused  by  geomagnetic  activity,  often  in  the  form 
of  geomagnetic  storms.  Ionospheric  storm  behavior  is  highly  variable,  because  of  the  several 
sources  of  enhancements  and  depletions.  General  morphological  studies  have  been  done 
although  individual  events  may  be  far  from  the  norm 

Generally,  at  midlatitudes  the  F i-layer  peak  density  undergoes  an  initial  enhance- 
ment (positive  phase),  followed  by  a sharp  drop  to  below  the  "quiet"  average  density.  The 
following  several  days  and  nights  may  remain  enhanced  or  depleted.  The  positive  phase  may 
be  explained  by  electrodynamic  lifting  (the  theory  of  equatorward  winds  blowing  the  ions 
up  the  tield  lines  has  lost  favor)  of  the  ions  into  lower-loss  regions,  causing  density  increases. 
The  negative  phase  may  be  either  a pushing  down  of  the  ions  into  higher-loss-rate  regions;  or, 
as  will  be  suggested  here,  the  propagating  of  turbulent  waves  out  of  the  polar  region  toward 
the  equator,  causing  increased  mixing  of  the  neutral  atmosphere,  which  in  turn  will  yield 
an  increased  loss  rate  at  all  altitudes. 

Davies  ( 1974)  suggested  that  there  is  a heat  source  on  the  auroral  oval  which  has  an 
apparent  localized  motion  that  is  greater  than  the  propagating  velocity  of  the  disturbance. 

This  gives  rise  to  a wake,  or  "shock  front,”  that  propagates  equatorward  and  results  in  a 
change  in  the  chemical  concentration  at  F-layer  heights.  The  neutral  concentrations  will 
return  slowly  to  undisturbed  conditions  over  a period  > 24  hours,  which  constitutes  the  “nega- 
tive phase”  of  an  ionospheric  storm. 

In  this  paper  we  shall  show  calculations  of  the  shock  front  shape  for  various  condi- 
tions an  outcome  of  this  model  is  that  a localized  source  area  (on  the  oval)  is  responsible 
for  the  portion  of  the  shock  front  at  an  observer’s  location.  In  other  words,  an  energetic 
event  along  southern  Greenland  may  affect  the  Eastern  United  States,  leaving  the  remaining 
portion  of  the  world  unchanged.  These  localized  storms  are  very  real  (Park,  1974),  and, 
using  this  present  theory,  we  may  learn  to  forecast  them. 


THE  MODEL 

The  Davies  model  considers  a point  heat  source  located  on  the  auroral  oval  (on  the 
noontime  magnetospheric  cusp)  that  moves  with  the  rotation  of  the  earth.  The  disturbance, 
here  presumed  to  be  gravity  waves,  propagates  radially  outward  at  a velocity  of  2°-5°  per 
hour.  Nighttime  joule  heating  would  have  similar  effects,  although  it  would  be  an  extended 
source. 

The  Feldstein  oval  is  placed  at  78°  (noon)  during  quiet  conditions,  and  Davies  uses  a 
2° -latitude  change  equatorward  per  unit  change  in  the  planetary  magnetic  index  K(V  Hence, 
the  source  point  is  at  a latitude  </>s: 

*s=78°-2Kp  . 


Considering  only  the  propagation  towards  lower  latitudes,  the  storm  front  will 
appear  as  a "wake”  similar  to  that  made  by  a speedboat.  As  it  moves  through  an  atmos- 
pheric region  (at  the  E-layer  input  altitude),  a cataclysmic  turbulent  mixing  takes  place, 
which  may  move  the  turbulent-mixing  boundary  (usually  at  an  altitude  of  100-120  km)  to 
heights  > 200  km.  Numerical  models  of  the  ionosphere  (Ruster  and  Dudeney,  1975)  indi- 
cate that  height  increases  in  this  boundary  level  result  in  large  decreases  in  F-region  electron 
density,  caused  by  changes  in  the  O/Nt  ratio.  Because  of  the  “shock  front"  moving  through, 
these  changes  can  occur  rapidly.  The  return  to  normal  quiet  conditions  will  be  at  the  slow 
rate  suggested  by  Rishbeth  and  Hanson  (1974).*  See  figure  Al. 


Figure  Al.  Geometry  of  the  shock  front  propagation 
calculation. 

0Q  is  the  source  latitude,  0 is  the  observer  latitude,  ^ is  the  angle  of  propagation 
from  the  normal. 

0 is  the  angle  swept  out  by  the  source  in  the  time  the  disturbance  propagates  a 
distance  8,  so  <M>0  is  the  longitude  difference  of  source  and  observation  point.  Notice  that 
it  is  latitude  dependent. 

The  propagation  path  length  8 is  described  by 

5 = a • <t>Q  , 

where  a = d8/d0  = cos  hence  ^ is  constant; notice  a = v/15,  where  v is  the  propagation 
velocity  (in  deg/h) 

Using  spherical  trigonometry, 

cos  0 = cos  0Q  cos  6 - sin  0Q  sin  6 cos  ( <p0  - <p) 


and 


•Rishbeth,  H,  and  Hanson,  WB,  “A  comment  on  plasma  ‘pileup'  in  the  F-region,"  J Atinos  and  Terr  Phys 
36,  703,  1974 


_.  sin  (s-0)  sin  (s-0o) 
<t>o  $ - ,an  sin  (s-8)  sin  (s) 


whore  s = (0  + 0Q  + 5 )/2  . 


Since  we  are  interested  in  the  longitude  (0o-0)  that  affects  a given  latitude  with  a 
source  at  (0O,  0O),  we  are  given  0,  0Q,  a,  cos  i p = 1 -a-.  Using  the  spherical  trigonometry 

law  of  cosines 


cos  0 = cos  0O  J\  - sin  “6  - sin  0Q  sin  6 cos  i p 


cos  0 + sin  0Q  sin  6 cos  4j  = cos  0Q  >/l  —sin-  . 

Squaring  and  regrouping: 

sin-8  jcos-  0Q  + sin-  0O  cos-  i|/  j + 
sin  8 ( 2 cos0  sin  0O  cos  \p  | + 

cos-  0 - cos-  0Q  = A sin “6  + B sin  6 + C = 0 

Solving  for  sin-  8 and  using  the  Law  of  Sines  ( sin  S/sin  ( <t>0~<P ) = sin  0/cos  i//).  we  get 

, , . v _ sin  8 / 2" 

sm  f^-^-^Vl-a 


sin  8 = 


B+v/B--4  AC 


The  two  roots  indicate  the  fact  that  there  is  poleward  propagation  as  well  as  equator- 
ward.  so  we  choose  the  positive  root,  as  we  are  interested  in  midlatitudes.  The  negative  root 
corresponds  to  over-the-pole  propagation  to  midlatitude. 

By  evaluating  8,  and  <t>0  - <t>  as  a function  of  latitude,  one  obtains  the  shock  front 
curve,  and  also  the  path  along  which  the  shock  front  “appears"  to  propagate  from  the  source 
equatorward.  Figure  A2  shows  the  fronts  and  paths  for  the  propagation  velocities  2°.  3°.  4°. 
5°/h.  The  spiral  structure  of  the  fronts  is  evident,  but  notice  that  an  energy  input  for  a short 
period  will  cause  only  a portion  of  this  spiral  to  develop  (fig  A3).  That  portion  will  propagate 
equatorward  along  the  propagation  paths,  so  by  setting  the  current  velocity  path  on  the 
observer’s  location,  one  can  determine  the  high-latitude  region  of  importance.  Greenland  will 
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affect  the  Atlantic  seacoast,  while  the  western  (ireat  Lakes  develop  shock  fronts  that  sweep 
across  the  California-Hawaii  Pacific  area. 

DISCUSSION 

The  outcome  of  this  model  (the  Davies  model)  is  that  a mixing  front  should  propagate 
equatorward  from  the  auroral  zone,  with  a rapid  depletion  of  the  F-region  (caused  by  greatly 
enhanced  loss  rates)  being  observed.  The  mixing  effectively  raises  the  turbopause,  or  lowers 
the  boundary  of  diffusive  equilibrium,  which  increases  the  molecular  nitrogen  concentration 
relative  to  atomic  oxygen;  this  increases  the  effective  loss  rate  at  F->  peak  altitudes,  followed 
by  a decrease  in  the  peak  electron  density.  Therefore,  if  the  model  is  valid,  a latitude  depend- 
ence in  the  “sudden  depletion”  time  should  be  observed,  and  should  fit  the  predicted  times 
of  arrival. 

Perhaps  more  important  is  the  conclusion  that  a specified  longitude  along  the  auroral 
oval  is  responsible  for  the  portion  of  the  front  passing  over  any  other  latitude,  longitude  posi- 
tion, and  that  energy  deposited  in  this  region  is  expected  to  affect  only  a localized  path 
(fig  A3). 

The  8 March  1970  ionospheric  storm  has  had  a complete  data  compilation  from  US 
stations  measuring  the  ionospheric  total  electron  content  (TKC)  (Klobuchar  and  others,  1971*). 
which  roughly  relates  to  F->  peak  electron  density.  One  set  of  the  observing  stations  lies  along 
a line  corresponding  to  our  postulated  "propagation  path,"  and  so  the  propagation  time  can 
be  checked  and  calculated.  Another  set  of  the  stations  lies  approximately  along  a constant 
geomagnetic  latitude  line,  so  an  estimate  of  the  longitude  extent  of  the  front  can  be  made. 
Figure  A4  shows  the  deviations  from  the  mean  of  the  electron  content  for  the  storm  day  for 
both  latitude  and  longitude  ordering.  Notice  that  along  the  propagation  path  the  stations  do 
show  increasing  delays  with  decreasing  latitude.  Figure  A5  compares  these  onset  times  with 
model  predictions  (for  various  propagating  velocities),  with  a velocity  ~3'/2°/h  fitting  the  data 
closest.  The  longitude  chain  does  indicate  a strong  dependence  on  longitude,  with  the  eastern 
stations  having  large  depletions  at  the  predicted  time,  and  towards  the  west  there  are  pro- 
gressively smaller  effects,  until  Hawaii  has  no  discernible  decrease.  This  is  expected  because 
the  Eastern  section  was  in  late  morning  conditions,  while  Hawaii  and  Stanford  were  presunrise; 
if  the  energy  is  dumped  in  the  noontime  cusp  for  a short  period,  the  effects  should  not  extend 
into  the  sunrise  longitudes. 

For  this  single  event  (the  most  thoroughly  documented)  the  model  provides  an  excel- 
lent description;  that  it  will  work  this  well  in  all  cases  is  unlikely,  but  more  extensive  testing  is 
necessary.  A major  problem  exists  in  studying  ionospheric  storms.  There  are  many  competing 
and  cooperating  processes  that  act  in  varying  magnitudes  to  make  each  event  individual.  Un- 
fortunately these  very  differences  impede  any  statistical  evaluation  of  ionospheric  storms  in 
general.  For  example,  if  the  “positive  phase"  is  extraordinarily  large,  the  depleting  effects  of 
one  front  (as  well  as  downward  motion)  may  leave  a net  increase  in  the  Fi  peak,  and  so, 
according  to  some  researchers,  no  negative  phase  occurred.  The  model  described  here  predicts 
sudden  depletions  in  the  F-region,  with  no  regard  to  the  initial  or  final  value  of  F-region 
density.  Any  other  approach  would  be  misleading. 


•Klobuchar,  JA.  and  others,  “Ionospheric  storm  of  March  8,  1970,"  J Geophys  Res,  76,  6202.  1971 
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Figure  A4.  Deviations  of  TEC  from  mean  for 
the  8 March  1970  event  showing  both  latitude 
and  longitude  ordering. 


Figure  AS.  Onset  time  (shock  front  passage) versus 
latitude,  with  data  of  8 March  1970  event  for 
comparison. 


SUMMARY 

The  "shock  front”  turbulent  mixing  model  initially  developed  by  Davies  and  investi- 
gated in  this  report  can  explain  some  of  the  phenomena  exhibited  in  ionospheric  storms.  In 
general,  ionospheric  storms  are  found  often  to  be  localized  in  impact,  with  sudden  decreases  in 
the  Fi-region  density  extending  far  down  in  latitude.  The  localized  storm  may  be  due  to  short 
and  localized  energetic  heating  in  the  auroral  zone,  which  causes  a disturbance  propagating 
equatorward  in  a narrow  “cone.”  The  disturbance  appears  to  propagate  with  a velocity  of 
~4°/h,  which  is  much  slower  than  the  Alfven  velocity  (the  velocity  at  which  any  plasma  dis- 
turbance should  propagate);  therefore,  the  disturbance  is  expected  to  be  propagating  through 
the  neutral  atmosphere,  perhaps  in  the  form  of  gravity  wave  trains.  That  gravity  waves  are 
caused  by  auroral  activity  is  well  documented  (Akasofu,  1968).  In  this  report  we  have  con- 
sidered the  possibility  of  developing  a shock-front-like  disturbance;  this  is  expected  because 
the  earth’s  rotational  motion  causes  the  auroral  oval  (and  polar  cusp)  to  move  at  supersonic 
velocities  relative  to  the  earth's  atmosphere.  The  model  does  describe  accurately  some  of  the 
events  seen  in  the  8 March  1970  event  documented  by  Klobuchar  and  others  (1971). 
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